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SUMMARY 


Fires  burning  simultaneously  over  hundreds  of  square  kilometers 
could  be  one  result  of  a  nuclear  weapon  explosion.  The  strong 
buoyancy  field  of  such  large  area  fires  induces  high-velocity  fire 
winds  that  turn  upward  in  the  burning  region  and  vertically  transport 
a  large  quantity  of  water  vapor.  In  this  report,  we  numerically  model 
the  rise  of  moisture-laden,  free-convection  columns  and  examine  the 
development  of  clouds  as  a  function  of  relative  humidity,  fire  size, 
and  burning  or  heat  release  rate.  The  rise  is  controlled  principally 
by  the  fire  heat  release  and  the  atmospheric  stratification.  In  most 
of  our  simulations,  enough  moisture  is  lofted  to  form  large  cumulus 
clouds  early  in  the  plume  development. 
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PREFACE 


This  work  was  sponsored  by  the  Defense  Nuclear  Agency  under 
contract  DNA  001 -85-C-0089.  The  work  was  monitored  by  Dr.  Michael  J 


CONVERSION  TABLE 


Conversion  factors  for  U.S.  Customary  to  metric  (SI)  units  of  measurement 


h  BY 

^  TO  OFT 

TA  r.tT  ^ _ _ _ 

RY^ _ 

_  DIVIDE 

angstrom 

1.000  000  X  E  -10 

meters  (m) 

atmosphere  (normal) 

1  013  25  X  E  *2 

kilo  pascal  (kPa) 

bar 

1.  000  000  XE  42 

kilo  pascal  (kPa) 

bam 

1.000  000  X  E  -28 

2  2 
meter  (m  ) 

British  thermal  unit  (thermochemical) 

1.  054  350  SE.3 

joule  (J) 

calorie  (thermochemical) 

4.  184  000 

Joule  (J) 

cal  (thermochemical) /cm 

4.  184  000  X  E  -2 

mega  joule/m2  (MJ/m2) 

curie 

3  700  000  X  E  4  1 

•giga  becquerel  (GBq) 

degree  (angle) 

1.  745  329  X  E  -2 

radian  (rad) 

degree  Fahrenheit 

=  Iff  *  459.  67)/l.  8 

degree  kelvln  (K) 

electron  volt 

1.  602  19  X  E  -19 

joule  (J) 

erg 

1.  000  000  X  E  -7 

joule  (J) 

erg /second 

1.  000  000  X  E  -7 

watt  (W) 

foot 

3.  048  000  X  E  -1 

meter  (mj 

f  oot  -pound  -f  o  rce 

1.  355  818 

joule  (J) 

gallon  (C  S  liquid) 

3  785  412  X  E  -3 

3  3 

meter  (m  ) 

inch 

2.  540  000  X  E  -2 

meter  (m) 

jerk 

1  000  000  X  E  *9 

joule  (J) 

joule /kilogram  (J /kg)  (radiation  dose 
absorbed) 

1.000  000 

Gray  (Gy) 

kilotons 

4  183 

terajoules 

kip  (1000  lbf) 

4.448  222  X  E  *3 

newton  (N) 

ktp/inch2  (ksi) 

6  894  757  XE  43 

kilo  pascal  (kPa) 

ktap 

l  000  000  X  E  +  2 

2 

newton  -second  /m 
(N-s/m2) 

micron 

1  000  000  X  E  -6 

meter  (m) 

mil 

2  540  000  X  E  -5 

meter  (m) 

mile  (international) 

1. 609  344  X  E  43 

meter  (m) 

ounce 

2  834  952  X  E  -2 

kilogram  (kg) 

pound -force  (lbs  avoirdupois) 

4.  448  222 

newton  (N) 

pound -force  inch 

1.  129  848  X  E  -1 

newton-meter  (N-  m) 

pound  -force /inch 

1  751  268  X  E  42 

newton /meter  (N/m) 

pound  -fo  rce  /Toot2 

4.  788  026  X  E  -2 

kilo  pascal  (kPa) 

pound  -force /inch"  (pai) 

6.  894  757 

kilo  pascal  (kPa) 

pound-mass  (tbm  avoirdupois) 

4  535  924  X  E  -1 

kilogram  (kg) 

pound -mass -foot2  (moment  of  inertia) 

4  214  011  X  E  -2 

kilogram-meter 
(kg  •  m2) 

pound -m ass  /Toot 

l  601  846  X  E  4  1 

kilogram  /meter 
(kg/m3) 

rad  (radiation  dote  absorbed) 

1  000  000  X  E  -2 

••Gray  (Gy) 

roentgen 

2  579  760  X  E  -4 

coulomb /kilogram 
(C/kg)  I 

shake 

1  000  000  X  E  -8 

aecond  (a)  ^ 

slug 

1  459  390  X  E  + 1 

kilogram  (Hg) 

torr  (mm  Hg.  0*  C) 

i — — — —  —  -J 

1  333  22  X  E  -1 

kilo  paacal  (kPa)  j 

•The  becquercl  (Bq)  is  the  SI  unit  of  radioactivity;  1  Bq  *  1  event/s. 
••The  Gray  (Gy )  is  the  S!  um(  of  absorbed  radiation. 


TABLE  OF  CONTENTS 


Section  Page 

SUMMARY  .  iii 

PREFACE  .  iv 

CONVERSION  TABLE  .  v 

LIST  OF  ILLUSTRATIONS  .  vii 

1  INTRODUCTION  .  1 

2  MODEL  DESCRIPTION  .  4 

Numerical  model  .  4 

Governing  equations  .  5 

Initial  and  boundary  conditions  . 6 

Heat  and  smoke  source  terms  .  7 

Water  phase  transitions  .  8 

3  RESULTS  .  10 

Input  parameters  .  10 

Baseline  water  cloud  .  11 

Baseline  smoke  cloud  .  14 

Comparison  between  dry  and  moist  atmosphere  .  15 

Effect  of  fire  area  and  heating  rate  .  17 

Fire  winds  .  17 

Effect  of  grid  resolution  .  20 

Individual  cases  .  22 

it  DISCUSSION  .  62 

5  LIST  OF  REFERENCES  .  64 


Appendices 

A  SMOKE  PRODUCTION  .  67 

B  WATER  PROCESSES  .  69 

C  WATER  CONTENT  DUE  TO  COMBUSTION  .  73 


vi 


! 


■■Ij'l.l-tJ'M  tJ’IJ  U  ‘■■  ‘■I  ' 


LIST  OF  ILLUSTRATIONS 


Figure 


Moisture  cloud  simulation  for  0.5=kW/m3  heating  rate 
and  7-km  radius  fire  . 


Smoke  cloud  simulation  for  0.5-kW/m3  heating  rate  and 
7~km  radius  fire  . 


3  Comparison  of  moist  and  dry  smoke  clouds  for  0.5-kW/m^ 
heating  rate  and  7-km  radius  fire  at  60  min  . 


Comparison  of  smoke  clouds  obtained  at  60  min  for  constant 
heating  and  variable  radius  and  for  constant  radius  and 
variable  heating  . 

Radial  velocity  at  8-km  radius  for  three  heating  rates 
for  7-km  radius  fire  at  10  min  intervals  . 


Smoke  clouds  at  45  min  for  fire  7  km  in  radius  with 
0.50  kW/m3  volume  heating  rate  . 


7  Moisture  cloud,  streamlines,  and  neutral  buoyancy  isotherms 
for  case  1 — high  heating  rate  (1.00  kW/m3),  large  radius 

( 1  0  km)  . 

8  Moisture  cloud,  streamlines,  and  neutral  buoyancy  iso¬ 

therms  for  case  2 — medium  heating  rate  (0.50  kW/m^),  large 
radius  ( 1 0  km)  . 

9  Moisture  cloud,  streamlines,  and  neutral  buoyancy  isotherms 
for  case  3 — low  heating  rate  (0.25  kW/m^),  large  radius 

(10  km)  . 

10  Moisture  cloud,  streamlines,  and  neutral  buoyancy  isotherms 

for  case  4 — high  heating  rate  (1.00  kW/m^),  medium  radius 
(7  km)  . 

11  Moisture  cloud,  streamlines,  and  neutral  buoyancy  isotherms 

for  case  5 — medium  heating  rate  (0.50  kW/m3),  medium  radius 
(7  km)  . 

12  Moisture  cloud,  streamlines,  and  neutral  buoyancy  isotherms 
for  case  6  —  low  heating  rate  (0.25  kW/m^),  medium  radius 

(7  km)  . 

13  Moisture  cloud,  streamlines,  and  neutral  buoyancy  isotherms 

for  case  ''--medium  hearing  ratD  (0.50  kW/m3),  small  radius 
(5  km)  . 


.v.v.v 


Witt 


1 

HI 


L  *.1  i'll1! 


Streamlines  and  neutral  buoyancy  isotherms  for  case  8-- 
dry  atmosphere  . 

Moisture  cloud,  streamlines,  and  neutral  buoyancy  iso¬ 
therms  for  case  9 — low- resolution  grid  . 


Smoke  cloud  and  velocity  vectors  for  case  1--high  heating 
rate  (1.00  kW/m^),  large  radius  (10  km)  . 


Smoke  cloud  and  velocity  vectors  for  case  2 — medium 
heating  rate  (0.50  kW/nP),  large  radius  (10  km) 


Smoke  cloud  and  velocity  vectors  for  case  3 — low  heating 
rate  (0.25  kW/m^),  large  radius  (10  km)  . . 


Smoke  cloud  and  velocity  vectors  for  case  4--high  heating 
rate  (1.00  kW/m3),  medium  radius  (7  km)  . 


Smoke  cloud  and  velocity  vectors  for  case  5 — medium 
heating  rate  (0.50  kW/m^),  medium  radius  (7  km)  .... 


Smoke  cloud  and  velocity  vectors  for  case  6--low  heating 
rate  (0.25  kW/m^),  medium  radius  (7  km)  . . 


Smoke  cloud  and  velocity  vectors  for  case  7--medium 
heating  rate  (0.50  kW/rtP),  small  radius  (5  km)  . 


Smoke  cloud  and  velocity  vectors  for  case  8--dry 
atmosphere  . . . 


Smoke  cloud  and  velocity  vectors  for  case  9--low-r esolution 


/.W/V  W  V  V  -J ir* 


SECTION  1 
INTRODUCTION 


A  nuclear  burst  over  an  urban  area  is  likely  to  cause  a  large 
number  of  fires.  The  initial  distribution  of  ignitions  is  related  to 
the  thermal  radiation  output  of  the  bomb  and  the  disruptive  effects 
(secondary  ignitions)  of  the  blast  wave.  The  precise  distribution  of 
fire  starts  depends  on  the  city  structure,  weapon  yield,  and  height  of 
burst.  It  is  reasonable  to  assume  a  frequency  of  ignitions  decreasing 
with  distance  from  the  burst.  For  a  1-MT  weapon,  ignitions  are  likely 
as  far  as  10  km  from  the  burst;  for  smaller  yields,  the  range 
decreases. 

Initially,  many  thousands  of  fire  starts  would  be  distributed 
across  an  urban  area.  Not  all  buildings  would  be  involved  simul¬ 
taneously,  but  it  is  likely  that  all  would  eventually  burn.  That  was 
suggested  by  the  area  fires  that  occurred  at  Hamburg,  Dresden, 
Hiroshima,  etc.  [Bond,  1946J.  During  the  course  of  an  area  fire, 
which  may  last  several  hours,  some  structures  would  be  actively  burn¬ 
ing,  some  would  be  in  preliminary  stages  of  involvement,  and  others 
would  have  already  completely  burned  out.  The  spread  responsible  for 
the  continual  burning  (or  heat  release)  of  the  city  is  limited  to  the 
area  of  fire  starts.  That  is,  a  large  area  is  not  burned  by  a  front 
sweeping  across  a  city,  but  rather,  by  an  evolving  complex  of  fires. 

Energy  would  be  released  over  a  large  area  rather  than  in  a 
narrow  fire  line,  generating  a  broad-based  convective  flow  or  plume. 
This  flow  would  be  supported  by  high-velocity  fire  winds  induced  by 
the  fire  complex  [Larson  and  Small,  1982;  Small,  Larson,  and  Brode, 
1984;  Small  and  Larson,  1984/85].  The  interactions  are 
straightforward.  Buoyancy  created  by  heat  release  generates  pressure 
gradients,  which  in  turn  induce  fire  winds.  Those  turn  upward  to  form 
the  convective  flow.  Radiative  cooling  rapidly  reduces  (over  ~10’  m) 
the  temperature  above  the  heat  release  region.  The  net  effect  is  a 
high-momentum  convective  flow  with  only  a  moderate  temperature  excess 


(several  degrees)  or  density  deficit  from  the  ambient.  For  axisym- 
metric  flows,  the  differential  is  largest  near  the  centerline. 

The  initial  momentum  and  temperature  excess  of  the  plume  depend 
directly  on  the  fire  size  and  heat  release.  These  determine,  in  part, 
the  plume  rise  and  deposition  of  smoke  in  the  atmosphere.  Entrainment 
of  ambient  air  during  the  rise  plays  a  minor  role  in  diffusing  the 
plume  momentum  and  buoyancy,  since  the  mixing  takes  place  at  the  plume 
edge,  far  from  the  largest  concentrations  of  buoyancy.  Consequently, 
atmospheric  structure  (lapse  rate,  inversion  height,  and  upper-level 
winds)  is  the  primary  factor  in  determining  the  ultimate  plume  rise. 
Fires  from  a  nuclear  burst  can  create  plume  motions  that  extend  to  the 
tropopause  [Hassig  and  Rosenblatt,  1983;  Small,  Remetch,  and  Brode, 
1984,  1985;  Penner,  Haselman,  and  Edwards,  1986;  Bacon,  Sarma,  and 
Proctor,  1986;  Cotton,  1985;  and  Tripoli  and  Kang,  1987]. 

The  high-velocity  fire  winds  induced  by  a  large  area  [0(100)  km2] 
city  fire  would  entrain,  at  low  altitude,  significant  amounts  of  water 
vapor.  Except  for  extremely  dry  atmospheres,  the  entrained  moisture 
is  roughly  an  order  of  magnitude  greater  than  the  water  produced  by 
combustion.  In  the  plume,  large  cloud  structures  evolve  from  conden¬ 
sation  of  the  entrained  moisture.  The  release  of  latent  heat  drives 
the  plume  higher;  evaporation  decreases  the  buoyancy.  Cloud  formation 
by  the  fire-generated  convective  motion  occurs  rapidly.  Such  cloud 
development  is  common  in  naturally  occurring  wildland  fires,  and  rain 
often  results.  Large  area  city  fires  would  produce  even  larger  and 
more  intense  clouds. 

Condensation  and  cloud  formation  not  only  influence  the  plume 
motion,  but  also  enable  the  scavenging  of  smoke  by  water  droplets,  and 
to  some  extent,  ice.  Early  cloud  formation  may  account  for  some 
scavenging  of  smoke  with  particularly  high  soot  content,  produced  soon 
after  the  fires  start.  Precipitation  may  limit  the  amount  of  smoke 
lofted  to  high  altitude.  Black  rain  was  noted  at  Hiroshima  [Committee 
for  the  Compilation  of  Materials  on  Damage  Caused  by  the  Atomic  Bombs 
in  Hiroshima  and  Nagasaki,  1981],  implying  that  fire-generated 
precipitation  removed  at  least  some  of  the  smoke. 
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The  rate  and  extent  of  cloud  formation  depends  on  the  low-level 
entrainment  of  humid,  ambient  air  and  it3  distribution  in  the  plume. 
Both  the  fire  winds  and  plume  structure  depend  directly  on  the 
buoyancy  field  or  details  of  the  fire.  Numerical  simulation  of  cloud 
formation  thus  requires  a  grid  fine  enough  to  properly  resolve  the 
fire-generated  flow,  as  well  as  the  larger  scale  atmospheric  response. 
In  this  report,  we  consider  two-dimensional  ( ax i symmetric)  approxima¬ 
tions  for  fire-generated  flows  in  a  stratified,  moist  atmosphere.  We 
have  determined  in  detail  the  amount  of  water  lofted,  plume  motion, 
cloud  development,  and  smoke  transport  for  several  fire  sizes  and 
burning  rates. 
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SECTION  2 
MODEL  DESCRIPTION 

NUMERICAL  MODEL. 

Our  numerical  simulation  employs  a  Lagrangian-Eulerian  algorithm 
to  model  atmospheric  convection  driven  by  a  large  area  fire,  and  is 
applicable  to  compressible,  two-dimensional  fluid  flows  over  a  wide 
range  of  flow  speeds  [Hirt,  Amsden,  and  Cook,  1974].  Compressibility 
is  important  near  the  fire  because  of  the  high  temperatures.  The 
Lagrangian-Eulerian  technique  is  a  compromise  between  fully  explicit, 
compressible  calculations  and  partially  implicit  calculations  as¬ 
sociated  with  the  anelastic  approximation  [Ogura  and  Phillips,  1962] 
for  deep  convection.  The  Lagrangian  phase  advects  the  mesh  vertices 
with  the  fluid,  assuming  no  flux  across  cell  boundaries.  In  the 
Eulerian  phase,  the  vertices  are  restored  to  their  original  positions 
and  fluxes  into  the  cells  are  accounted  for.  Upwind,  or  donor  cell, 
differencing  is  used  in  the  Eulerian  phase;  centered  differences  in 
the  Lagrangian  phase.  The  technique  allows  high  resolution  solutions 
using  variable  zoning  of  the  comoressible  flow  near  the  fire  and 
within  the  plume  and,  at  the  same  time,  it  permits  time  steps  much 
larger  than  possible  with  purely  explicit  schemes. 

The  algorithms  used  in  this  study  include  many  enhancements  to 
the  original  Lagrangian-Eulerian  scheme  used  by  Hirt,  Amsden,  and  Cook 
[1974].  Pressure  is  calculated  implicitly  at  the  end  of  each  time 
step.  Temperature  T  is  determined  from  the  energy  equation,  which 
includes  a  radiation  term  proportional  to  T1*  and,  therefore,  the 
temperature  is  also  found  implicitly.  We  have  also  modified  the 
original  program  to  include  latent  heating  or  cooling  due  to  water 
vapor  phase  transitions,  thermal  diffusion,  volumetric  heating  in  the 
fire  region,  and  the  advection  of  water  and  smoke. 

Variable  zoning  is  used.  Fine  resolution  is  achieved  near  the 
lower  boundary,  the  centerline,  and  the  fire  region.  A  fixed  number 
of  grid  points  (19  radial,  including  the  axis;  4  vertical,  including 
the  ground)  is  usually  used  within  the  fire  volume  independent  of  fire 


radius.  The  spacing  increases  by  a  factor  of  approximately  1.1  for 
cells  outside  the  fire  region.  The  axial  grid  spacing  within  the  fire 
volume  is  33  m»  increasing  to  about  1  km  at  the  tropopause.  The 
computational  domain  extends  vertically  to  30  km  and  radially  to 
70  km. 


GOVERNING  EQUATIONS. 

The  simulations  are  based  on  the  compressible,  nonhydrostatic 
Navier-Stokes  equations  with  velocity  components  (u,v),  in  the  radial 
and  axial  directions  (r,y): 


3u  ^  3u  3u 
3t  U  3r  V  3y 


3v  3v  3v 
3t  +  U  3r  +  V  3y 


i  +  k  /  3_  1  3ru  +  3  u 
p  3r  l  3r  r  3r  ^ 2 


1  3d  „  /  1  3  3v 

p  3y  g  +  KM  l  r  3r  r  3r 


where  t  is  time,  p  is  pressure,  is  eddy  viscosity,  and  g  is 
gravitational  acceleration.  The  ideal  gas  equation  of  state  is  used. 
Although  velocities  are  subsonic,  0(1)  changes  in  temperature  and 
density  occur  as  a  result  of  the  large  heat  release  within  the  fire 
volume.  The  full  continuity  equation  is  thus  required: 


Ifi  +  1  (rpu) 
3t  r  3r  vrpu' 


The  energy  equation  is 


3c  T  3c  T  3c  T 
v  v  v 

3t  +  u  3r  +  V  3y  ' 


-  BT  (i  p- 
\r  3r 


where  q^.,  q^,  and  q^  are  heat  sink/ source  terms  for  heating  by  the 
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fire,  radiative  cooling,  and  latent  heating  due  to  changes  in  water 
phase,  respectively.  The  gas  constant  R  and  specific  heat  volume  cy 
are  mass-weighted  properties: 


p  "  pd  +  pv  +  P1  +  pi  ’ 


R  = 


p  R  +  p  R  +  p  R,  +  p.  R. 
d  c  v  v  11  i  l 


P  c  +  p  c  p,  c  ,  +  p.  c  , 
c  vd  v  vv  1  vl  i  vi 


(4) 


The  total  mixture  density  is  p  and  the  subscripts  d,  v,  1,  and  i  refer 
to  dry  air,  water  vapor,  liquid  water,  and  ice,  respectively. 

INITIAL  AND  BOUNDARY  CONDITIONS. 

Normal  velocities  are  zero  at  all  boundaries.  A  no-slip  boundary 
condition  is  applied  at  the  ground  and  free-slip  conditions  are  used 
at  the  axis  of  symmetry,  upper  boundary,  and  outer  radius.  It  is 
assumed  that  the  upper  boundary  is  isothermal;  the  lower  boundary  and 
centerline  are  perfect  insulators.  The  atmosphere  is  initially  quies¬ 
cent,  and  the  temperature  is  prescribed  by  the  U.S.  Standard  Atmos¬ 
phere  [1962],  The  initial  lapse  rate  is  stable  everywhere  and  con¬ 
stant  (-6.45  x  10“3°c/m)  fran  the  ground  to  an  altitude  of  10  km  in 
the  troposphere,  isothermal  ( —5 6°C )  in  the  tropopause  from  10  to 
20  km,  and  constant  (1.06  x  10~3°C/m)  from  20  to  30  km  in  the 
stratosphere.  The  stable  stratospheric  layer  at  the  top  and  the  large 
radial  extent  of  the  computation  domain  isolate  the  flow  from  the  top 
and  outer  boundaries.  An  initial  vertical  distribution  of  relative 
humidity  [Manabe  and  Wetherald,  1967]  is  assumed: 


't  the  ground,  4>0  is  77  percent  and  d0  is  the  ambient  surface  pres- 


HEAT  AND  SMOKE  SOURCE  TERMS. 

The  three  energy  source/sink  terms  (fire  heat  release,  radiative 

cooling,  and  latent  heat  release)  control  the  system  dynamics.  The 

fire  is  modeled  as  a  100  m-high,  volume  heat  source.  The  heating  rate 

increases  linearly  in  time  from  qf  =  0  to  a  constant  value  at  t  = 

900  s  and  thereafter.  Radius  and  heating  rate  are  varied  to  test  the 

effects  of  fire  size  and  intensity.  Initially,  all  properties  are 

horizontally  homogeneous,  except  within  the  fire  volume  where  eddy 

viscosity  and  eddy  thermal  conductivity  are  set  an  order  of  magnitude 

3  2 

larger  than  outside  the  fire  volume  (K  T  -  1 0  m  /s,  as  opposed  to 


K  _  =  10  m  /s).  The  larger  values  account  for  fire-generated  tur- 
M ,  T 

bulence  and  effective  surface  roughness  associated  with  urban 
features.  Smoke  generation  and  transport  are  described  by: 


♦  ;  h  <ruS)  *  »7  <vs) 


1  .66 


— 6  • 
X  10  q 


where  S  is  the  smoke  concentration,  and  S  is  the  source  function  based 

on  a  3  Percent  smoke  emission  rate  and  heat  release  rate  q„  (see 

2  1 

Appendix  A).  The  combustion  of  2  g/cm  /h  of  fuel  is  assumed  to 
2 

release  100  kW/m  .  Heat  release  increases  as  fuel  loading  increases; 
the  total  smoke  mass  input  to  the  atmosphere  depends  on  the  burnable 
fuel  loading.  During  the  time  of  significant  heat  release  (roughly 
0.5  to  4  h) ,  smoke  transport  is  primarily  determined  by  large-scale 
features  of  the  velocity  field,  and  therefore,  is  modeled  as  an  advec- 
tive  process  [see  Eq.  (6)].  Diffusion  becomes  a  significant  transport 
process  only  at  much  later  times. 

Our  use  of  a  volume  heat  source  to  model  the  distribution  of 
fires  results  in  an  average  temperature,  in  the  region  of  heat 
release,  considerably  lower  than  that  of  the  individual  fires.  Ac- 
cordingly,  the  radiative  loss  qp  is  modeled  by: 


Jn 


g 


(7) 


Here,  a  is  the  ratio  of  actual  burning  area  to  the  total  area,  H  is 
the  height  of  the  fire  region,  T0  is  the  far-field  ambient  tempera¬ 
ture,  and  o  is  the  Stef an-Boltzmann  constant.  In  general,  structures 
occupy  less  than  30  percent  of  the  land  in  dense  urban  areas.  Not  all 
the  area  burns  at  one  time,  nor  does  all  of  the  area  contain  acces¬ 
sible  combustible  material.  Thus  a  <  1 .  The  ratio  B  of  actual  fire 
temperature  to  average  region  temperature  is  roughly  3.  Temperature 
decays  rapidly  above  the  fire  as  a  result  of  radiative  cooling. 

WATER  PHASE  TRANSITIONS. 

A  significant  amount  of  heat  is  released  by  water  phase  transi¬ 
tions  in  the  plume  well  above  the  fire;  the  effect  of  latent  heat 
release  becomes  increasingly  important  for  lower  fire  heating  rates. 

At  each  time  step  in  the  calculation,  moist  air  is  first  transported 
by  the  flow  without  phase  change.  The  saturation  vapor  density  is 
then  compared  to  the  new  water  vapor  density  at  each  grid  point.  A 
phase  transition  then  occurs,  bringing  the  vapor  density  as  close  as 
possible  to  saturation  by  means  of  evaporation,  condensation,  melting, 
freezing,  or  sublimation,  depending  on  temperature  (see  Appendix  B) . 
For  temperatures  warmer  than  0°C,  the  process  requires  that,  first, 
all  ice  melts,  and  then,  either  liquid  water  evaporates  or  water  vapor 
condenses,  contingent  on  whether  the  air  is  subsaturated  or  super¬ 
saturated  with  respect  to  liquid  water.  At  temperatures  lower  than 
-15°C,  liquid  water  first  freezes,  and  then  ice  sublimates  or  vapor 
accretes,  depending  on  whether  the  air  is  subsaturated  or  super¬ 
saturated  with  respect  to  ice. 

For  tne  temperature  range  0°  to  -15°C,  liquid  water  and  ice  are 
allowed  to  coexist  to  simulate  the  supercooled  water  condition  ob¬ 
served  in  naturally  occurring  clouds.  A  deficit  in  water  vapor  den¬ 
sity  relative  to  saturation  over  ice  is  first  corrected  by  evaporating 
water,  since  the  saturation  vapor  pressure  over  water  is  higher  than 
that  over  ice,  and  then  any  further  subsaturation  is  made  up  for  by 
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sublimating  ice.  Alternatively,  an  excess  water  vapor  density  rela¬ 
tive  to  saturation  over  water  is  corrected  by  the  accretion  of  vapor 
to  ice.  If  the  water  vapor  pressure  is  subsaturated  with  respect  to 
liquid  water  but  saturated  with  respect  to  ice,  no  change  in  vapor 
density  takes  place.  Thus,  liquid  water  may  be  advected  into  regions 
where  the  temperature  is  in  the  range  0°  to  -15°  C  but  does  not  form 
in  these  regions.  This  is  consistent  with  the  fact  that  the  vapor 
pressure  over  water  is  greater  than  that  over  ice  within  this  tempera¬ 
ture  range.  All  temperature  changes  due  to  latent  heating  or  cooling 
are  assumed  to  occur  at  constant  pressure  and  are  incorporated  in  the 
energy  equation  at  the  next  time  step  [6t  is,  at  most,  0(1 0~ 1 )  s]. 

Upon  condensation  (or  freezing),  one-third  (or  one-half)  of  the 
condensate  is  removed  frcm  the  atmosphere  to  simulate  precipitation. 
Details  of  the  precipitation  process,  particularly  scavenging  of  smoke 
particles,  are  not  included  in  this  preliminary  study.  This  modeling 
of  precipitation,  although  it  is  rudimentary,  does  reproduce  the 
irreversible  process  found  in  naturally  occurring  clouds  of  moist 
adiabatic  ascent  followed  by  dry  adiabatic  descent,  warming  the  cloud 
region. 

Combustion  of  the  hydrogen  component  of  the  fuel  in  the  presence 
of  atmospheric  oxygen  produces  a  mass  of  water  approximately  half  the 
fuel  mass.  This  water  vapor  is  also  drawn  up  into  the  plume,  where  it 
condenses  and  releases  latent  heat.  However,  the  mass  of  water  vapor 
produced  by  combustion  is  generally  at  least  an  order  of  magnitude 
lower  than  that  present  in  the  atmosphere  or  entrained  by  the  fire 
winds.  Except  for  a  very  dry  atmosphere,  it  has  a  negligible  effect 
on  the  plume  dynamics  (see  Appendix  C). 


SECTION  3 
RESULTS 


INPUT  PARAMETERS. 

The  test  matrix  of  input  conditions  covered  a  range  of  fire 
sizes,  fire  intensities,  and  atmospheric  moisture  conditions  (see 
Table  1).  Area  fires  of  radii  5,  7,  and  10  km  were  considered.  They 
correspond  either  to  cities  completely  involved  in  fire  but  of  dif¬ 
ferent  sizes,  or  to  fires  started  by  weapons  of  different  yields. 
Heating  by  the  fire  does  not  take  place  directly  at  the  surface,  but 
is  distributed  across  a  volume  of  finite  depth  as  a  result  of  combus¬ 
tion  of  fuel  vapors  within  the  volume,  and  conduction  and  radiation  of 
heat  away  from  localized  fire  regions.  A  1 00-m  depth  was  assumed. 
Combustion  of  2  g/cm2  (1  to  2g/cm2  is  a  typical  fuel  loading  for 
residential  areas)  in  1  h  releases  100  kW/m2,  or  1  kW/m^  throughout 
the  heat  release  volume.  One-  to  four-hour  burn  times  are  reasonable; 
higher  fuel  loadings  and  longer  burn  times  are  plausible.  A  7-km 
radius,  0.5-kW/m3  fire  would  be  typical  for  most  city  sizes,  fuel 
loadings,  and  expected  burning  rates. 


Table  1.  Test  matrix  and  summary  of  input  conditions. 


Heating  Rate 

Fire  Radi 

Case 

(kW/m^ ) 

(km) 

/  1 

1.00 

1  0 

1  2 

0.50 

10 

1  3 

0.25 

1  0 

Moist  /  4 

1 .00 

7 

1  5 

0.50 

7 

1  6 

0.25 

7 

V  7 

0.50 

5 

Dry  8 

0.50 

7 

Grid 


9 


0.50 


7 
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Simulations  were  also  performed  for  fire  intensities  of  0.25, 
0.50,  and  1.00  kW/m^.  These  values  represent  either  different  fuel 
loadings  or  burn  times.  With  the  exception  of  one  simulation  in  a  dry 
atmosphere,  all  cases  had  identical  atmospheric  temperature  and  mois¬ 
ture  structure.  In  addition,  one  case  tested  the  effect  of  grid 
resolution.  Results  include  smoke  concentrations,  distributions  of 
liquid  water  and  ice,  and  velocity,  pressure,  temperature,  and  density 
fields . 

BASELINE  WATER  CLOUD. 

Many  features  of  plume  evolution  are  similar  for  all  the  fire 
intensities  and  sizes  tested.  Figures  1  and  2  show  the  development  of 
the  moisture  and  smoke  clouds  for  a  baseline  case  (7~km  radius  fire, 
0.5~kW/m3  heat  release  rate).  Initially,  there  is  a  warming  of  the 
lowest  2  km  of  the  atmosphere  above  the  fire.  Air  drawn  radially 
inward  toward  the  fire  produces  an  upward  mass  flux.  That,  in  turn, 
produces  subsidence,  slightly  warming  the  entire  troposphere  at  larger 
radii.  The  inflowing  air  initially  turns  up  sharply  near  the  outer 
perimeter  of  the  fire  and  forms  a  small  cloud  (Fig.  la).  The  cloud 
base  forms  at  2  km  but  rises  slightly  as  the  fire  (heat  release)  warms 
the  lower  atmosphere.  Gradually,  the  converging  radial  inflow  estab¬ 
lishes  a  strong  updraft  at  the  centerline  (Fig.  1c).  Later,  as  the 
buoyant  air  gains  momentum,  it  is  lofted  considerably  above  its  level 
of  neutral  buoyancy  and  becomes  colder  and  more  dense  than  the  ambient 
air  surrounding  the  plume  (Figs.  le-f).  This  cold  air  diverges 
radially  at  the  too  of  the  plume,  falls  downward,  again  overshoots 
equilibrium,  and  becomes  warmer  than  the  surrounding  air.  The  process 
may  continue  several  times  at  successively  larger  radii.  Thus,  the 
cloud  consists  of  a  warm  core,  sheathed  by  alternating  cold  and  warm 
regions.  The  core  is  cloud  free  near  the  base  because  the  air  is 
sufficiently  warm  to  be  subsaturated.  Similar  cloud-free  vaults  are 
commonly  observed  in  large  thunderstorms.  The  low  pressure  caused  by 
the  heat  of  the  fire  produces  subsidence  and  warming  at  high 
altitudes.  Local  thickening  of  the  stratosphere  may  result. 
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Figure  2.  Smoke  cloud  simulation  for  0.5-kW/m-^  heating  rate  and  7-km 
radius  fire. 


Early  (less  than  20  min)  in  the  evolution  of  the  flow  field,  a 
small  vortex  develops  at  the  fire  edge.  Vertical  velocities  in  the 
plume  core  are  influenced  by  the  inward  translation  of  that  vortex. 

The  sharp  thermodynamic  gradients  at  the  fire  edge  cause  a  periodic 
reforming  of  the  vortex  near  the  fire  perimeter,  reducing  the  flow  to 
the  centerline.  As  a  result,  the  vertical  velocity  of  the  plume  core 
oscillates.  Between  25  and  40  min,  the  plume  top  remains  relatively 
stationary  at  about  14.5  km  as  a  new,  secondary  circulation  forms  near 
the  perimeter  of  the  fire  (Figs.  Ic-e).  When  this  pattern  breaks 
down,  allowing  increased  inflow  to  the  centerline,  the  plume  top  rises 
to  17  km  (not  shown  in  Figs.  1  or  2)  at  55  min.  Later,  as  the  outer 
secondary  circulation  reemerges,  plume  core  velocities  decrease  and 
the  top  drops  to  about  15.5  km.  The  centerline  oscillations  cause 
large  amplitude  gravity  waves  that  propagate  radially  outward 
(Fig.  Ic-f). 

BASELINE  SMOKE  CLOUD. 

Figure  2  shows  the  smoke  plume  evolution  for  the  same  radius  and 
heating  rate  as  the  moisture  cloud  sequence  in  Fig.  1.  The  smoke  and 
moisture  clouds  coincide,  except  near  the  ground  and  centerline,  where 
the  rising  warm  air  is  subsaturated,  and  near  the  lower,  outer  edge, 
where  smoke  detrains  from  the  moisture  cloud  (see  Figs.  If  and  2f ) . 

The  top  of  the  smoke  cloud  approximately  coincides  with  the  top  of  the 
ice  cloud  at  the  top  of  the  dome  of  cold  air  capping  the  warm  plume 
core.  The  cold  dome  represents  the  maximum  vertical  extent  reached  by 
air  from  the  surface,  before  its  density  excess  causes  it  to  fall  back 
to  lower  altitudes.  There  is  a  well-defined  boundary  between  the 
smoke  cloud  and  the  clear  air  outside  it.  The  strongest  concentra¬ 
tions  of  smoke  lie  near  the  source  region  and  the  plume  core.  The 
axisymmetric  geometry  of  the  simulation  must  be  taken  into  account 
when  considering  smoke  distribution.  In  Fig.  2f ,  over  70  percent  of 
the  total  smoke  mass  lies  outside  the  4-km  radius,  despite  the  higher 
concentrations  for  smaller  radii. 

Large-scale  mixing  processes  spread  the  smoke  out  from  the  plume 
core  at  higher  altitudes.  The  vortex  above  the  outer  edge  of  the  fire 
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appears  to  fold  pockets  of  clear  ambient  air  into  the  smoke  cloud 
(Fig.  2e).  The  vortex  also  seems  to  produce  a  downward  rolling  motion 
near  the  outer  edge.  Such  features  are  consistent  with  observations 
of  actual  fire  plumes.  Except  for  the  largest,  most  intense  fires, 
smoke  does  not  appear  to  be  injected  into  the  stratosphere.  The 
maximum  altitude  reached  by  the  smoke  from  a  0.5-kW/m^,  7-km  fire  was 
1 7  km  at  55  min.  Simulations  of  more  intense  fires  (1  kW/m3)  result 
in  correspondingly  higher  plume  rises  (discussed  later);  lower  heat 
release  rates  (0.25  kW/m3)  result  in  lower  plume  heights.  The  outcome 
is  more  or  less  independent  of  fire  area,  although  a  threshold  size 
presumably  exists  [Small  and  Larson,  1984/5].  Comparable  plume  heights 
have  been  observed  in  nature  [Manins,  1985].  The  radial  spreading  of 
the  smoke  is  independent  of  fire  radius  or  heat  release  rate,  and  is 
roughly  confined  to  a  region  less  than  twice  the  radius  of  the  fire  at 
60  min.  Upper-level  winds  would  presumably  transport  smoke  downwind, 
or  diffusion  would  spread  smoke  radially,  at  later  times.  However, 
the  large  radial  spread  predicted  by  Penner,  Haselman,  and  Edwards 
[1986]  in  quiescent  atmospheres  is  not  reproduced  in  these  simula¬ 
tions  . 

COMPARISON  BETWEEN  DRY  AND  MOIST  ATMOSPHERE. 

Comparison  of  plume  growth  in  a  moist  atmosphere  to  that  in  a  dry 
atmosphere  shows  that  smoke  deposition  i3  markedly  affected  by  ambient 
moisture.  Figure  3  shows  the  difference  in  plume  growth  after  60  min. 
The  heat  release  rate  is  the  same  as  the  baseline  case  shown  in 
Figs.  1  and  2  (0.5  kW/m^,  7-km  radius).  When  plume  motions  reach  the 
level  of  condensation  and  latent  heat  is  released,  buoyancy 
increases.  Comparison  of  the  maximum  centerline  velocities  shows  the 
effect  of  moisture.  The  velocity  for  the  moist  case  (Fig.  3a)  is 
117  m/s  and  for  the  dry  case  (Fig.  3b)  it  is  81  m/s.  Condensation 
produces  high-velocity  motions  early  in  the  development  of  the  plume, 
carrying  large  concentrations  of  smoke  to  high  altitudes.  The  smoke 
spreads  out  radially  initially  by  advection  and  later  by  diffusion. 

The  maximum  vertical  and  radial  extent  of  the  smoke  clouds  are  about 
the  same,  but  the  moist  case  deposits  more  smoke  in  the  upper  atmos¬ 
phere. 
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EFFECT  OF  FIRE  AREA  AND  HEATING  RATE. 


fires  of  constant  radii  but  three  different  heat  release  rates  are 


increased  heat  release  rate  for  fires  of  the  same  radius.  Clearly, 


classical  plume  theories  are  unsuitable  for  modeling  large-scale 


a  finite  area. 


FIRE  WINDS. 


Figure  5  shows  the  fire  wind  velocity  at  8  km  for  the  7-km 
baseline  radius  (1  km  outside  the  fire  volume)  and  the  three  heating 


rates.  Maximum  inflow  velocities  occur  at  early  times  and  near  the 


jjj 


The  simulations  clearly  show  condensation  early  in  the  plume 
evolution.  Figure  1  shows  initial  condensation  at  15  min  (Fig.  la) 
and  shows  ice  developing  within  25  min  (Fig.  1c).  Ice  formation  is 
generally  indicative  of  the  onset  of  precipitation.  Although  not 
explicitly  modeled,  precipitation  early  in  the  plume  development 
suggests  that  seme  of  the  dense  sooty  smoke  released  by  the  fire  in 
its  early  stages  may  be  scavenged.  That  appears  to  be  consistent  with 
observations  of  rain  generated  by  wildland  fires  and  of  black  rain 
generated  by  several  fires  during  World  War  II  [Committee  for  the 
Compilation  of  Materials  on  Damage  Caused  by  the  Atomic  Bombs  in 
Hiroshima  and  Nagasaki,  1981].  The  efficiency  of  smoke  particles, 
however,  as  cloud  condensation  nuclei,  is  not  well  known,  and  only  a 
small  fraction  may  actually  be  scavenged  [Pittock,  et  al.,  1986]. 


Clouds  produced  by  fires  of  three  different  radii  but  equal  heat 
release  rates  are  shown  in  Figs.  4a  through  4c,  and  clouds  produced  by 


shown  in  Figs.  4d  through  4f.  Figure  4  shows  that  a  concentrated  heat 
source  deposits  smoke  much  higher  in  the  atmosphere  than  does  a  source 
of  equal  volume-integrated  magnitude  that  is  spread  over  a  larger 
surface  area.  Figure  4a  through  4c  shows  that  plume  height  at  60  min 
is  roughly  the  same  (it  actually  decreases  slightly  from  17.5  to 
15.5  km)  as  the  fire  radius  increases,  even  though  total  heat  release 
increases  by  a  factor  of  4.  In  contrast,  Figs.  4d  through  4f  show 
that  plume  height  increases  (from  9  to  23  km)  in  direct  response  to 


fires — they  use  a  point  source,  rather  than  distributing  energy  across 
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Figure  4.  Comparison  of  smoke  clouds  obtained  at  60  min  for  constant 
heating  and  variable  radius  and  for  constant  radius  and 
variable  heating. 
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ground.  Later,  the  inflow  wind  speed  is  lower  but  occupies  a  deeper 
layer,  so  the  depth- integrated  inflow  remains  nearly  constant  after 
25  min.  There  is  outflow  aloft  associated  with  the  return  circula¬ 
tions  shown  in  Fig.  1.  The  outflow  is  most  intense  at  early  times 
when  the  azimuthally  aligned  vortex  lies  just  over  the  heating  volume 
circumference.  The  mean  fire  wind  velocities  near  the  ground  are 
surprisingly  low  in  comparison  to  observations  of  fire  winds  as¬ 
sociated  with  large  urban  fires  [Committee  for  the  Compilation  of 
Materials  on  Damage  Caused  by  the  Atomic  Bombs  in  Hiroshima  and 
Nagasaki,  1981 ;  Small  and  Larson,  1984].  The  maximum  is  20  m/s  but 
the  average  is  approximately  10  m/s.  Observations  of  fire  winds  are 
highly  subjective  and  reported  wind  speeds  are  not  usually  based  on 
instrumented  measurements.  Winds  may  be  channeled  by  streets  and 
structures — such  winds  may  not  be  representative  of  the  mean  flow. 
Moreover,  the  model  is  two-dimensional  and  does  not  simulate  channeled 
winds  nor  those  associated  with  local  concentrations  of  vorticity. 

EFFECT  OF  GRID  RESOLUTION . 

The  effect  on  the  baseline  case  of  a  reduction  of  grid  resolution 
was  tested  by  changing  the  radial  spacing  from  an  average  spacing  of 
about  270  m  near  the  centerline  (increasing  by  15  percent  for  succes¬ 
sive  grid  points  in  the  positive  radial  direction)  to  a  constant 
1  km.  The  vertical  spacing  was  changed  from  10C  m  near  the  ground, 
(increasing  by  20  percent  for  successive  grid  points  above  the  ground) 
to  1000  m  (increasing  by  20  percent  for  successive  grid  points  above 
the  ground).  The  lower  resolution  case  is  comparable  to  that  used  in 
three-dimensional  simulations.  Figure  6  shows  the  effect  of  reduced 
resolution  on  the  baseline  smoke  cloud  (Fig.  6a).  Clearly,  smoke  is 
not  lofted  as  high  in  the  low- resolution  case  as  it  is  in  the  high- 
resolution  case  for  identical  heating  and  fire  size.  There  appears  to 
be  a  greater  concentration  of  smoke  at  low  altitude  for  the  low- 
resolution  case.  The  reasons  for  such  dramatic  differences  are  un¬ 
clear,  but  the  higher  resolution  is  expected  to  more  accurately  simu¬ 
late  the  thermodynamics  of  smoke-cloud  evolution.  Lower  resolution, 
two-  and  three-dimensional  results  [Penner,  Haselman,  and  Edwards, 
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1986;  and  Tripoli  and  Kang,  1987]  in  fact,  do  sometimes  show  greater 
horizontal  smoke  spreading  and  lower  lofting. 

INDIVIDUAL  CASES. 

In  addition  to  the  general  comparisons  of  cases  with  differing 
heating  rates,  radii,  and  moisture  content  in  this  section,  we  high¬ 
light  some  features  of  the  individual  moist-atmosphere  simulations. 
Plots  of  water  and  smoke  clouds  are  shown  in  Figs.  7  through  15  and  16 
through  24,  respectively. 

Case  1:  High  Heating  Rate  (1.00  kW/m^),  Large  Radius  (10  km). 

A  strong  circulation  initially  occurs  at  the  fire  edge  and  lasts 
20  min.  During  that  period,  there  is  little  penetration  radially  to 
the  centerline  and  little  vertical  motion  to  the  condensation  level. 
Maximum  lofting  into  the  tropopause  occurs  at  35  min.  After  that 
time,  two  low-level,  azimuthally  aligned  vortices  appear  within  the 
troposphere  and  divert  the  flow  from  the  centerline.  The  fire  ini¬ 
tially  heats  the  lower  levels.  The  rising  air  overshoots  its  level  of 
equilibrium  and  becomes  cooler  than  the  surrounding  air.  Due  to 
adiabatic  compression,  air  in  the  downward  leg  of  the  circulation 
becomes  warmer  than  ambient.  The  fire  continually  heats  the  lower 
levels,  particularly  near  the  centerline.  Shear  stress  appears  to 
carry  air  aloft  (despite  its  negative  buoyancy)  at  the  plume  edge.  A 
weak  secondary  circulation  within  the  upper  tropopause  and  lower 
stratosphere  carries  air  downward  near  the  centerline  and  is  as¬ 
sociated  with  the  low  pressure  caused  by  the  fire  heat  near  the  sur¬ 
face  . 

Case  2:  Medium  Heating  Rate  (0.50  kW/m^),  Large  Radius  (10  km). 

Streamline  gradients  near  the  ground  indicate  large  mass  flux  and 
high  wind  into  the  fire  region.  Cloud  formation  begins  at  20  min. 

The  cloud  base  is  elevated  near  the  centerline,  forming  a  cloud-free 
vault  after  30  min.  Again,  the  formation  of  the  double  azimuthally 
oriented  vortex  configuration  corresponds  to  a  decrease  in  radial 
inflow  to  the  centerline.  The  temperature  perturbation  from  ambient 
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Figure  7.  Moisture  cloud,  streamlines,  and  neutral  buoyancy  isotherms 
for  case  1  - -h i gh  heating  rate  (1.00  kW/m^),  large  radius 
km) 
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Figure  7.  Moisture  cloud,  streamlines,  and  neutral  buoyancy  isotherms 
for  case  1  - -h i gh  heating  rate  (1.00  kW/rn^),  large  radius 
(10  km)  (Concluded). 
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Figure  8. 


Moisture  cloud,  streamlines,  and  neutral  buoyancy  isotherms 
for  case  2--medium  heating  rate  (0.50  kW/rn^),  large  radius 
(10  km). 
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Figure  8. 


Moisture  cloud,  streamlines,  and  neutral  buoyancy  isotherms 
for  case  2--medium  heating  rate  (0.50  kW/rn-^),  large  radius 
(10  km)  (Concluded). 
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Figure  10.  Moisture  cloud,  streamlines,  and  neutral  buoyancy  isotherms 
for  case  4--high  heating  rate  (1.00  kW/m^),  medium  radius 
(7  km)  (Concluded). 
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Figure  14.  Streamlines  and  neutral  buoyancy  isotherms  for  case  8--dry 
atmosphere  (Concluded). 
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Figure  15.  Moisture  cloud,  streamlines,  and  neutral  buoyancy  isotherms 
for  case  9--low-resolution  grid. 
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Figure  16. 


Smoke  cloud  and  velocity  vectors  for  case  1 --high  heating 
rate  (1.00  kW/m3),  large  radius  (10  km)  (Concluded). 
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Figure  17.  Smoke  cloud  and  velocity  vectors  for  case  2--medium 

heating  rate  (0.50  kW/nw),  large  radius  (10  km)  (Concluded) 
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Figure  18,  Smoke  cloud  and  velocity  vectors  for  case  3--low  heating 
rate  (0.25  kW/m3),  large  radius  (10  km). 
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Figure  19<  Smoke  cloud  and  velocity  vectors  for  case  4 - -h i gh  heating 
rate  (1.00  kW/m-3),  medium  radius  (7  km). 
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Figure  19.  Smoke  cloud  and  velocity  vectors  for  case  4--high  heating 
rate  (1.00  kW/m^),  medium  radius  (7  km)  (Concluded). 
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Figure  20.  Smoke  cloud  and  velocity  vectors  for  case  5--medium 

heating  rate  (0.50  kW/rrP),  medium  radius  (7  km)  (Concluded). 
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Figure  21.  Smoke  cloud  and  velocity  vectors  for  case  6--low  heating 
rate  (0.25  kW/m-*),  medium  radius  (7  km). 
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Figure  21.  Smoke  cloud  and  velocity  vectors  for  case  6--1ow  heating 
rate  (0.25  kW/rn^),  medium  radius  (7  km)  (Concluded). 
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Figure  22.  Smoke  cloud  and  velocity  vectors  for  case  7--medium 
heating  rate  (0.50  kW/nv5),  small  radius  (5  km). 
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Figure  22.  Smoke  cloud  and  velocity  vectors  for  case  7--medium 

heating  rate  (0.50  kW/rrw),  small  radius  (5  km)  (Concluded] 
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is  positive  (warm)  near  the  cloud  core,  negative  within  a  region 
capping  and  surrounding  the  warm  central  core,  and  then  positive  again 
in  the  region  capping  and  surrounding  the  aforementioned  cold  region. 
Apparently,  the  warm  air  rising  along  the  centerline  (the  first  warm 
region)  overshoots  its  equilibrium  level  because  of  its  momentum.  This 
air  continues  to  rise,  but  due  to  moist  adiabatic  expansion  it  becomes 
colder  than  its  surroundings  (the  first  cold  region).  It  then  falls 
downward,  but  once  again  overshoots  equilibrium  and  warms  due  to 
adiabatic  compression  (the  second  warm  region).  This  process  may  be 
repeated  at  successively  larger  radii  and  corresponds  to  gravity  waves 
propagating  radially  outward.  The  base  of  the  cloud  initially  con¬ 
sists  entirely  of  liquid  water  but  later  becomes  a  mixture  of  super¬ 
cooled  water  and  ice.  The  cloud  base  first  occurs  near  the  0°C  level 
(2.3  km)  but  later  rises  to  as  high  as  5  km  due  to  low-level  heating. 
At  later  times  (50  to  55  min)  an  observed  feature  of  large  smoke 
clouds  is  reproduced — the  folding  into  the  smoke  column  of  pockets  of 
clear  air  near  the  base  of  the  mushrocm  cap.  This  process  apparently 
enhances  the  mixing  of  clear  surrounding  air  with  the  smoke  by  bring¬ 
ing  clear  air  into  the  vicinity  of  greatest  turbulence  and  smoke 
concentration. 

Case  3:  Low  Heating  Rate  (0.25  kW/m3) ,  Large  Radius  (10  km). 

This  case  has  the  lowest  heat  release,  buoyancy  production,  and 
thus,  the  smallest  pressure  gradients.  The  fire  edge  vortex  is  ac¬ 
cordingly  sustained  longer  than  in  the  two  previous  cases,  and  the 
water  and  moisture  clouds  never  reach  significant  height  (9  km 
maximum).  For  this  "weak"  heating  the  flow  develops  slowly  and  the 
mass  flux  to  the  centerline  is  limited.  The  plume  and  cloud  also 
develop  slowly. 

Case  High  Heating  Rate  (1.00  kW/m^),  Medium  Radius  (7  km). 

Strong  updrafts  associated  with  the  concentrated  heat  source  near 
the  surface  carry  air  and  smoke  above  the  level  of  neutral  buoyancy 
within  10  min.  The  warm-core,  cloud-free  vault  is  visible  at  all 
times,  not  just  in  the  later  stages  of  cloud  development.  Pockets  of 
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clear  air  are  constantly  being  advected  into  the  cloud.  The  top  of 
the  smoke  cloud  reaches  22  km — the  lower  stratosphere.  The  inflow 
boundary  layer  is  unusually  thick.  The  large  radial  flux  of  air  into 
the  fire  region  is  consistent  with  the  lofting  of  air  to  high  al¬ 
titude. 

Case  5:  Medium  Heating  Rate  (0.50  kW/m^),  Medium  Radius  (7  km). 

This  is  the  baseline  case  previously  discussed.  The  complete 
time  sequence  is  given  in  Figs.  11  and  20. 

Case  6:  Low  Heating  Rate  (0.25  kW/m^),  Medium  Radius  (7  km). 

The  radial  extent  of  the  fire  is  too  great  to  sustain  a  large 
organized  circulation  to  the  centerline.  Soon  after  the  flow  reaches 
the  centerline,  the  vortex  at  the  fire  perimeter  reforms,  reducing 
flow  to  the  inner  radii.  The  cloud  top  seems  to  be  stationary  at 
1 2  km  by  55  min,  and  the  only  change  thereafter  is  the  radial  spread 
of  smoke  at  lower  altitudes.  Much  of  the  inflow  between  1  and  3  km 
appears  to  be  driven  by  viscous  entrainment  (which  "drags"  the  air 
inward  over  the  fire  region)  rather  than  by  pressure.  This  inflow  is 
at  such  a  high  altitude  that  it  is  not  heated  by  the  fire.  Thus,  it 
never  has  much  buoyancy  and  does  not  develop  sufficient  momentum  to  be 
lofted  to  higher  altitudes.  It  quickly  returns  to  its  equilibrium 
level.  This  is  clearly  evident  in  the  streamline  patterns  at  late 
times. 

Case  7:  Medium  Heating  Rate  (0.50  kW/m^),  Small  Radius  (5  km). 

The  characteristic  pattern  of  a  warm  core  cloud  surrounded  by 
alternating  regions  of  cold  and  warm  air  occurs  within  40  min.  The 
inflow  velocity  achieves  a  maximum  of  15  m/s  (54  km/h).  Between  55 
and  60  min,  a  cold  bubble  of  air  separates  from  the  cloud  at  about  the 
level  of  the  cloud  base.  At  the  same  time,  a  pocket  of  clear  air  is 
being  folded  into  the  core  region  of  the  cloud.  The  cold  bubble  is 
probably  a  result  of  the  evaporation  of  condensate  as  the  air  detrains 
from  the  cloud,  since  the  flow  is  nearly  horizontal  and  does  not 
appear  to  undergo  adiabatic  expansion.  This  process  is  reminiscent  of 
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smoke  detraining  from  the  base  of  moisture  clouds  formed  by  large 
wildland  fires:  The  downward  rolling  motion  near  the  cloud  base 
causes  the  moisture  component  to  evaporate  leaving  cold,  stable,  smoky 
air,  which  then  separates  fran  the  cloud.  Since  it  is  well  above  the 
region  of  inflow,  it  does  not  immediately  reenter  the  smoke  column. 


Case  8:  Dry  Atmosphere. 

The  complete  time  sequence  for  the  dry  atmosphere  case  is  shown 
in  Figs.  1 4  and  23. 

Case  9:  Low  Grid  Resolution. 

The  complete  time  sequence  for  the  low  grid  resolution  case  is 
shown  in  Figs.  15  and  2^. 
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SECTION  4 
DISCUSSION 


The  motions  generated  by  large  area  fires  are,  in  general,  three- 
dimensional;  the  burning  region  is  likely  to  be  asymmetric,  and  wind 
shears  can  obviously  influence  the  plume  rise  and  downwind  diffusion 
of  smoke.  Nevertheless,  plume  motions — even  in  strong  winds — are 
often  virtually  axisymmetric  over  most  of  the  rise  in  the  atmosphere. 
The  approximation  is  best  near  the  fire  or  region  of  strong  buoyancy 
and  least  appropriate  at  the  plume  top.  Although  there  is  some  loss 
in  generality,  the  higher  resolution  of  an  axisymmetric  simulation 
more  realistically  models  (within  the  burning  region  and  plume  volume) 
the  entrainment  of  low-level  moisture,  the  water-phase  transitions, 
and  the  smoke  transport  than  does  a  three-dimensional  simulation  of 
lower  resolution. 

We  use  a  volume  heat  source  rather  than  a  surface  heat  flux  to 
model  the  fire  energy  release.  The  volume  source  is  an  approximation 
of  the  concentrated  effects  of  actual  fire  complexes,  but  it  does 
provide  a  more  realistic  distribution  of  buoyancy  than  does  a  surface 
condition.  As  a  consequence,  more  accurate  representations  of  the 
pressure  field,  fire  winds,  and  low-level  moisture  convergence  into 
the  plume  are  obtained.  Our  simulations  show  that  for  a  range  of  fire 
sizes  and  intensities,  plume  development  is  primarily  determined  by 
fire  intensity  and  is  little  affected  by  fire  area  (for  radii  greater 
than  5  km) . 

We  have  shown  that  atmospheric  moisture  contributes  significantly 
to  plume  evolution,  and  that  early  scavenging  of  smoke  particles  by 
precipitation  is  likely  to  reduce  the  amount  injected  into  the  upper 
atmosphere.  Moreover,  the  confinement  of  the  plume  to  a  local  area 
suggests  greater  smoke  deposition  locally.  Except  for  the  largest 
fires,  the  tropopause  is  an  effective  barrier  against  the  injection  of 
smoke  into  the  stratosphere;  since  the  height  of  the  tropopause  varies 
with  season  and  latitude,  these  features  must  be  taken  into  account 
when  determining  plume  injection  heights. 
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Rain  can  influence  the  plume  dynamics.  Rain  either  can  be 
evaporated  and  thus  cool  the  air  or  fall  to  the  ground.  Cooling  the 
air  above  the  fire  reduces  its  buoyancy  and  inhibits  the  plume  rise. 
However,  latent  heat  is  stored  during  the  evaporation  process  and  its 
later  release  within  the  rising  plume  may  actually  enhance  the  plume 
rise.  The  fire-generated  cloud  extends  radially  beyond  the  fire 
region.  Precipitation  from  such  clouds  is  more  likely  to  evaporate 
before  reaching  the  ground  above  the  fire  than  outside  the  fire 
region.  Moisture  entrainment  by  the  fire  winds  to  the  burning  region 
could  thus  be  increased.  This  influences  the  energy-momentum  ex¬ 
changes  (perhaps  only  slightly).  In  some  cases,  fire  spread  is  pos¬ 
sible  beyond  the  initial  fire  start  area;  rain  could  prevent  or  at 
least  moderate  such  spread.  Smoke  scavenging  and  deposition  of 
radioactive  material  beyond  the  principal  weapon  effects  region  by 
fire-induced  (or  augmented)  rain  may  be  significant.  The  intei — 
relationship  between  thermodynamics,  combustion,  and  dynamics  is 
subtle  and  not  yet  completely  defined  for  large  area  fires;  further 
modeling  and  simulations  are  in  progress. 
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APPENDIX  A 
SMOKE  PRODUCTION 


The  smoke  production  rate  per  unit  area  S.  for  a  3  percent  smoke 

2  A 

emission  (assuming  a  2  g/cm  /h  burn  rate)  is 


=  burn  rate  x 


smoke  mass 
fuel  mass 


?  g  fuel  „  0.03  g  smoke 

2  1.00  g  fuel 

cm  h 


-0.166  8-5S2SS  . 

m  s 


Now  let  the  smoke  produced  at  the  surface  immediately  fill  a  100-m 
high  volume.  Then  the  smoke  release  rate  per  unit  volume  Sy  is 


-3  g  smoke 


S./100  m  =  1 .66  x  1 0 

A 


2  3 

Assuming  that  a  burn  rate  of  2  g/cm  /h  corresponds  to  a  1  kW/m 
volume  heating  rate  in  a  100-m  high  volume  and  that  burn  rate  is 
directly  proportional  to  a  heating  rate  qf,  the  volume  smoke  produc¬ 
tion  is  given  by 


Sv  =  S  =  1 .66  x  10 


smoke  t 
m^  s  kW/m^ 


S  -  1.66  x  10 


-6  kg  smoke 


m^  s  kW/m^ 


APPENDIX  B 
WATER  PROCESSES 


The  type  of  water  phase  transition  that  takes  place  is  determined 
by  the  temperature,  availability  of  condensate,  condensate  phase,  and 
vapor  density  relative  to  the  temperature-dependent  saturation  value. 
The  incremental  latent  heat  release  Aq,  during  a  single  time  step  At, 
is 

Aq  =  ~Llv  Apv  "  Lli  Api  =  cp  +  Ap/p  '  (12) 

where  L^y  is  latent  heat  of  vaporization  [0.^60^  x  10^  -  2.369  x  10^ 

(T  -  Tf)  J/kgJ,  is  latent  heat  of  fusion  (0.3336  x  10^  J/kg),  pv 
is  vapor  density,  is  ice  density,  cp  is  specific  heat  at  constant 
density,  Tf  is  freezing  point  temperature  (-15°C),  and  T  is 
temperature.  This  is  a  constant  pressure  process  and,  therefore,  Ap  = 
0.  The  change  in  water  vapor  density  occurring  during  At  is 

APy  =  min  (pvs  -  pv,  Pj  +  Pj)  ,  (13) 

where  Rv  is  the  gas  constant  for  water  vapor  (0.J46014  x  1 0 3  j/kg/K), 
pvs  is  saturation  vapor  pressure  or  610.78  exp  [6.8355  x  1C)3  (l/Tf  - 
1/T)  +  5.1^55  In  Tf/T],  pV3  is  saturation  vapor  density  (pvs  = 
Pyg/Ry/T),  and  p^  is  liquid  density.  Thus,  if  supersaturated  condi¬ 
tions  exist,  there  is  a  decrease  in  vapor  density  equal  to  - | pvs  - 
p y | .  If  subsaturated  conditions  exist,  the  water  vapor  deficit  | pvs  - 
pv|  is  compensated  for  by  a  conversion  of  condensate  to  vapor  until 
saturation  is  reached,  as  long  as  p^  +  p^  exceeds  | pvs  -  pv| .  If  | pV3 
-  Py |  exceeds  p^  +  p^,  the  increase  in  water  vapor  density  is  p^  +  p^. 

The  parameters  defining  the  amount  of  precipitation  are  set  so 
that  one-third  the  liquid  water  and  one-half  the  ice  formed  during  At 
are  removed  by  precipitation: 
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If  the  air  is  subsaturated  (Apv  £  0)  at  the  beginning  of  the  time  step 
then  there  is  no  precipitation,  so: 


a.  =  1  . 
l 


The  water  vapor  density  is  updated  and  the  preliminary  heat  release  is 


p  p  +  Ap 

V  V  rv 


q  <r-  -L  Ap  . 
lv  v 


where  q  is  heat  release  in  joules  per  cubic  meter.  The  decision¬ 
making  process  built  into  the  code  takes  account  of  the  temperature. 
If  the  temperature  exceeds  the  melting  point  Tm(=  0°C),  then  cooling 
takes  place  due  to  the  fact  that  all  ice  melts: 


q  q  -  l1.  p. . 


The  amount  of  liquid  water  increases  by  the  amount  of  melted  ice.  It 
decreases  by  evaporation  if  the  air  is  subsaturated  (Apv  >  0,  =  1) 

or  increases  by  the  amount  of  condensation  not  lost  to  precipitation 
if  the  air  is  supersaturated  (Apv  <  0,  =  2/3): 


P1  P1  +  pi  "  al  Apv’ 


and  finally 


Ki 
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p .  0  . 
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Due  to  the  scarcity  of  suitable  freezing  nuclei  in  the  atmosphere,  the 


freezing  point  Tf  is  taken  to  be  -15°C  rather  than  0°C.  Supercooled 
water  and  ice  coexist  within  the  temperature  range  Tm  £  T  Z  Tf.  A 
high  concentration  of  smoke  particles  may  alter  this  phenomenon  but 
barring  definitive  observations  to  the  contrary,  we  adhere  to  the 
stated  parameterization. 

If  Apv  £  pj_  Z  0  (subsaturated)  and  p^  2  0,  cooling  occurs  due  to 
sublimination  of  ice  and  evaporation  of  water: 


q  q  -  Lu(Apv  -  P]_ )  . 


The  term  Lf  f  (Apv  -  p^)  arises  from  the  water  vapor  deficit  being 
first  made  up  by  evaporation  of  liquid,  then  by  supercooled  water,  and 
finally  by  sublimation  of  ice.  This  is  commensurate  with  the  vapor 
pressure  over  water  being  less  than  that  over  ice  at  the  same 
temperature.  The  change  in  the  liquid  and  ice  components  is,  respec¬ 
tively. 


p.  p.  +  p.  -  Ap  , 
l  l  1  v 


P1  0  . 


If  Apv  S  pj_,  the  change  in  the  liquid  water  component  is 


°1  P1  '  al  &Pv  * 


Since  <*i  55  <*i  =  1  for  Apv  >  0,  the  process  used  for  the  supercooled 
water  temperature  range  also  holds  for  T  £  Tf. 
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APPENDIX  C 

WATER  CONTENT  DUE  TO  COMBUSTION 


The  amount  of  water  normally  present  in  the  atmosphere  far  ex¬ 
ceeds  that  produced  by  combustion.  Consider  the  combustion  of  F,  a 
typical  fuel  load  F  (F  =  10  kg/m2).  The  hydrogen  content  of  urban 
fuels  is  approximately  5  percent  by  mass,  slightly  higher  (-  6.5 
percent)  for  cellulose.  The  molecular  weight  of  hydrogen  is  1  and 
that  of  oxygen  is  16.  A  5  percent  hydrogen  content,  when  combined 
with  atmospheric  oxygen,  results  in  water  due  combustion  of  0.45  F. 
Assume  that  the  water  produced  at  the  surface  is  lofted  relatively 
slowly  at  vertical  velocity  v  (v  =  3  m/s),  so  as  to  maximize  water 
vapor  density,  and  that  it  fills  a  volume  directly  above  the  surface 
where  it  was  produced.  If  the  fuel  is  completely  consumed  in  time 
t  =  1  h,  then  the  water  vapor  density  pv  in  the  column  above  the  site 
of  combustion  is 

pv  =  0.45  F/v/t  =  4.2  x  10'4  ^  H20  .  (24) 

combustion  nr 

Now  consider  the  atmospheric  water  vapor  density  for  a  low  rela¬ 
tive  humidity  (20  percent) ,  so  as  to  minimize  the  atmospheric  water 
vapor  density.  The  relative  humidity  is  the  ratio  of  ambient  mixing 
ratio  w  to  the  saturation  mixing  ratio  u>s: 

mixing  ratio  =  —  =  —  =  0.2  ,  (25) 

U)  6 

S  S 

where  e  and  e3  are  the  ambient  and  saturation  vapor  pressures,  respe>.~ 


At  20°C , 


es(20°C)  =  23.373  *  IQ2  N/m2  , 


4.68  x  102  N/m2 


ambient 


RyT  ( 46o  J/kg/K)  (293  K) ’ 


=  3-5  x  1<)"3  kg/m3 


Thus,  even  under  extreme  conditions,  the  atmospheric  water  vapor 


density  exceeds  that  produced  by  combustion  by  an  order  of  magnitude. 


IB 

rsx. 


a 


$ 

85 


II 


‘  >1. 


1  Afm"  if | 


DISTRIBUTION  LIST 


DEPARTMENT  OF  DEFENSE 


DEPARTMENT  OF  THE  ARMY 


ARMED  FORCES  RADIOBIOLOGY  RSCH  INST 
ATTN:  VBOGO 

ASSISTANT  TO  THE  SECRETARY  OF  DEFENSE 

ATOMIC  ENERGY 

ATTN:  LTCOLL  MILLS 

DEFENSE  INTELLIGENCE  AGENCY 
ATTN:  DB-6E2  C  WIEHLE 
ATTN:  N BARON 
ATTN:  RTS-2B 
ATTN:  RTS-2 B 
ATTN:  WDB-4CR 

DEFENSE  NUCLEAR  AGENCY 
ATTN:  DFRA 
ATTN:  DFSPG  ULLRICH 
ATTN.  OPNS 
ATTN:  RAAE 

ATTN.  RAAE  K  SCHWARTZ 
ATTN.  RAAE  L  WITTWER 
ATTN.  RAEEG  BAKER 
ATTN:  RAEERWEBB 
ATTN:  RARPDAUTON 
ATTN:  SPAS  M  FRANKEL 
ATTN:  TDTD/C  CORSETTI 
2  CYS  ATTN:  TDTR 
4  CYS  ATTN.  TITL 


DEFENSE  TECHNICAL  INFORMATION  CENTER 
12CYS  ATTN:  DD 

DIRECTOR 

ATTN:  LTCOL  G  BETOURNE 
ATTN:  R  RUFFIN 

FIELD  COMMAND  DEFENSE  NUCLEAR  AGENCY 
ATTN:  FCTXE 
ATTN:  FTTDWSUMMA 

JOINT  STRAT  TGT  PUNNING  STAFF 
ATTN:  JKCS 

NATIONAL  DEFENSE  UNIVERSITY 
ATTN:  COLS  GARDINER 
ATTN:  G  FOSTER  MOB  CONCEPTS  DEV  CTR 
ATTN:  H  ALMOND 

OFFICE  OF  THE  SECRETARY  OF  DEFENSE 
ATTN:  COL  A  RAMSAY 


U  S  ARMY  ATMOSPHERIC  SCIENCES  UB 
ATTN:  RSUTHERUND 
ATTN:  SLCAS-AR-M  MR  RUBIO 

U  S  ARMY  CORPS  OF  ENGINEERS 

ATTN:  DAEN-RDM  R  GOMEZ 
ATTN:  DR  CHOROMOKOS  DAEN-RDM 

U  S  ARMY  CORPS  OF  ENGINEERS 
ATTN.  L  ZIEGLER 
ATTN:  R  BECKER 

U  S  ARMY  ENGR  WATERWAYS  EXPER  STATION 
ATTN:  L  LINK 

U  S  ARMY  MISSILE  INTELLIGENCE  AGENCY 
ATTN:  J  GAMBLE 

U  S  ARMY  NATICK  RSCH  DEV  &  ENGRG  CENTER 
ATTN:  HMEL-BISI 

U  S  ARMY  STRATEGIC  DEFENSE  COMMAND 
ATTN:  DR  J  LILLY 
ATTN:  GEDLIN 
ATTN:  J  VEENEMAN 
ATTN.  M  CAPPS 
ATTN:  R BRADSHAW 

DEPARTMENT  OF  THE  NAVY 

CNO  EXECUTIVE  PANEL 

ATTN.  CAP  L  BROOKS 

NAVAL  RESEARCH  UBORATORY 
ATTN:  RJEK 

NAVAL  SURFACE  WEAPONS  CENTER 
ATTN:  K-44S  MASTERS 

DEPARTMENT  OF  THE  AIR  FORCE 

AF/INYXC 

ATTN.  LTCOL  N  BARRY 

AIR  FORCE  GEOPHYSICS  UBORATORY 
ATTN:  D  CHISHOLM 
ATTN:  LS/R  MURPHY 
ATTN:  LSI/  H  GARDINER 
ATTN:  LYC/R  SANTA 
ATTN:  LYPHSMUENCH 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY/EN 
ATTN:  AFIT/ENP  MAJ  S  R  BERGGREN 


s? 

5$ 


Dist-l 


WWAiTi 


DNA-TR-87-176  (DL  CONTINUED) 


AIR  FORCE  OFFICE  OF  SCIENTIFIC  RSCH 
ATTN:  D  BALL 

AIR  FORCE  SPACE  DIVISION 

ATTN:  YNC  CAPT  K  O'BRYAN 

AIR  FORCE  TECHNICAL  APPLICATIONS  CTR 
ATTN:  J  MARSHALL 

AIR  FORCE  WEAPONS  LABORATORY 
ATTN:  CAPTLEONG 
ATTN:  JJANNI 
ATTN:  J  W  AUBREY,  NTED 

BALLISTIC  MISSILE  OFFICE 

ATTN:  LT  ROTHCHILD 

ATTN:  MYSP/CAP  TOMASZEWSKI 

DEPUTY  CHIEF  OF  STAFF/XOX 
ATTN:  AFXOX 

STRATEGIC  AIR  COMMAND/XPXF 
ATTN:  T  BAZZOLI 

DEPARTMENT  OF  ENERGY 

ARGONNE  NATIONAL  LABORATORY 
ATTN.  H  DRUCKFR 
ATTN:  M  WESLEY 

BROOKHAVEN  NATIONAL  LABORATORY 
ATTN.  B  MANOWITZ 
ATTN:  E  WEINSTOCK 

DEPARTMENT  OF  ENERGY 
ATTN:  I  NEDDOW 
ATTN:  T  HARRIS 

DESERT  RESEARCH  INSTITUTE 
ATTN:  JHALLETT 
ATTN:  J  HUDSON 

LAWRENCE  BERKELEY  NATIONAL  LAB 
ATTN.  H  ROSEN 

LAWRENCE  LIVERMORE  NATIONAL  LAB 
ATTN:  C  R  MOLENKAMP 
ATTN:  C  SHAPIRO 
ATTN:  F LUTHER 
ATTN:  GBING 
ATTN:  G  SIMONSON 
ATTN:  JPENNER 
ATTN:  J  POTTER 
ATTN:  L-10  A  GROSSMAN 
ATTN:  L-262  A  BROYLES 
ATTN:  L-262  J  KNOX 
ATTN:  L-442,  J  BACKOVSKY 
ATTN:  L-453  L  ANSPAUGH 
ATTN:  M  MACCRACKEN 
ATTN:  N  ALVAREZ 
ATTN:  R  MALONE 


ATTN:  RPERRET 
ATTN:  RPERRETT 
ATTN:  SGHAN 

LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  D  SAPPENFIELD 
ATTN:  DR.  D  CAGLIOSTRO 
ATTN:  E  J  CHAPYAK 
ATTN:  E  JONES 
ATTN:  E  SYMBALISTY 
ATTN:  G  GLATZMAIER/ESS 
ATTN:  GM  SMITH 
ATTN:  L  HAUER 
ATTN:  L  CLOUTMAN 
ATTN:  PHUGES 
ATTN:  R  MALONE 
ATTN:  T  YAMATTA 

OAK  RIDGE  NATIONAL  LABORATORY 
ATTN:  D  FIELDS 

SANDIA  NATIONAL  LABORATORIES 
ATTN:  A  L  JOHNSON 
ATTN:  BZAK 

ATTN:  D  DAHLGREN  6440 

ATTN:  D  FORDHAM 

ATTN:  D  WILLIAMS 

ATTN:  DIV  G-449  K  D  BERGERON 

ATTN:  LTROST 

ATTN:  MD  BENNETT 

ATTN:  ORG  332  R  C  BACKSTROM 

OTHER  GOVERNMENT 

CENTRAL  INTELLIGENCE  AGENCY 
ATTN:  AWARSHAWSKY 
ATTN:  7E47R  NELSON 

DEPARTMENT  OF  AGRICULTURE 
ATTN:  D  HAINES 

DEPARTMENT  OF  COMMERCE 
ATTN:  H  BAUM 
ATTN:  R  LEVINE 

DEPARTMENT  OF  TRANSPORTATION 
ATTN:  COL  M  ROESCH 

DIRECTOR,  FFASR 

ATTN:  C  CHANDLER 

ENVIRONMENTAL  PROTECTION  AGENCY 
ATTN.  R  COTHERN 
ATTN:  WE  FALLON 

FEDERAL  EMERGENCY  MANAGEMENT  AGENCY 
ATTN:  BW  BLANCHARD 
ATTN:  D  BENSON  NP-CP-MR 
ATTN:  DKYBAL 
ATTN:  HTOVEY 


DNA-TR-87-176  (DL  CONTINUED) 


ATTN:  J  POWERS 
ATTN:  J  RUMBARGER 
ATTN:  OFC  OF  RSCH/NP  H  TOVEY 
ATTN:  S  ALTMAN 

GENERAL  ACCOUNTING  OFFICE 
ATTN:  A  PIERCE 
ATTN:  PJBOLLEA 
ATTN:  V  BIELECKI 

NASA 

ATTN:  N  CRAYBILL 
ATTN:  WRCOFER 

NASA 

ATTN:  RHABERLE 
ATTN.  OTOON 
ATTN:  R  YOUNG 
ATTN:  T  ACKERMAN 

NATIONAL  BUREAU  OF  STANDARDS 
ATTN.  G  MULHOLLAND 
ATTN:  R  LEVINE 
ATTN:  RREHM 
ATTN:  R  SCHRACK 

NATIONAL  BUREAU  OF  STANDARDS 
ATTN.  H  BAUM 

NATIONAL  CENTER  ATMOSPHERIC  RESEARCH 
ATTN:  JKIEHL 
ATTN:  S  SCHNEIDER 
ATTN:  S  THOMPSON 
ATTN.  VRAMASWAMY 

NATIONAL  CLIMATE  PROGRAM  OFFICE 
ATTN:  AHECHT 
ATTN:  MYERG 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ATTN:  FFEHSENFELD 
ATTN:  JDELUISI 
ATTN:  R  PUESCHEL 
ATTN:  V  DERR 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN. 
ATTN:  B  HICKS 

NATIONAL  RESEARCH  COUNCIL 
ATTN:  KBEHR 
ATTN:  R  DEFRIES 

NATIONAL  SCIENCE  FOUNDATION 
ATTN:  EBIERLY 
ATTN:  HVIRJI 
ATTN:  LHAMATY 
ATTN:  R  SINCLAIR 
ATTN:  R  TAYLOR 

NUCLEAR  REGULATORY  COMMISSION 
ATTN:  R  ALEXANDER 


OFFICE  OF  SCIENCE  AND  TECH  POLICY 
ATTN:  BHEALY 
ATTN:  COLS  WYMAN 

OFFICE  OF  TECHNOLOGY  ASSESSMENT 
ATTN:  R  WILLIAMSON 

U  S  ARMS  CONTROL  &  DISARMAMENT  AGCY 
ATTN:  B  DOENGES  NWC-DPA 
ATTN:  G  PITMAN 
ATTN:  H  SCHAEFFER 
ATTN:  R  GODESKY 
ATTN:  R  O'CONNELL  NWC-DPA 

II  S  DEPARTMENT  OF  STATE 
ATTN:  ACORTE 
ATTN:  C  CLEMENT, (OES.STS) 

ATTN:  S  CLEARY 

ATTN:  T  VREBALOVICH  OES.STS 

US  GEOLOGICAL  SURVEY 
ATTN;  R  DECKER 

U  S  GEOLOGICAL  SURVEY 

ATTN:  E  SHOEMAKER 

U  S  HOUSE  OF  REPRESENTATIVES 
A1TN:  C  BAYER 

ATTN:  COMMITT  ON  SCI  &  TECH  J  DUGAN 

U  S  HOUSE  OF  REPRESENTATIVES 
ATTN:  J  FREIWALD 
ATTN:  MHERBST 

US  DEPARTMENT  AGRICULTURE 
ATTN:  DWARD 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AERO-CHEM  RESEARCH  LABS,  INC 
ATTN:  DB  OLSON 

AERODYNE  RESEARCH,  INC 
ATTN:  C  KOLB 
ATTN:  J  LURIE 

AEROJET  ELECTRO- SYSTEMS  CO 
ATTN:  AFYMAT 
ATTN:  S  HAMILTON 
ATTN:  R  PAN 

AEROSPACE  CORP 

ATTN:  CRICE 
ATTN;  LR  MARTIN 

AEROSPACE  CORPORATION 
ATTN:  G  LIGHT 

ALLEN  RESEARCH  CORP 
ATTN.  R  ALLEN 


DNA-TR-87-176  (DL  CONTINUED) 


AMERICAN  ASSN  ADVANCEMENT  OF  SCIENCE 
ATTN:  DM  BURNS 

ANALYTIC  SERVICES,  INC  (ANSER) 

ATTN:  RBROFFT 

APPLIED  RESEARCH  CORP 
ATTN.  AENDAL 

AT&T  DEFENSIVE  SYSTEMS  STUDIES 
ATTN:  RJANOW 

ATMOSPHERIC  AND  ENVIRONMENTAL  RES 
ATTN:  NSZE 

AUDIO  INTELLIGENCE  DEVICES  INC 
ATTN:  H  BAUM 

AVCO  CORPORATION 

ATTN:  G  GRANT,  DEPT  MGR 

BALL  AEROSPACE  SYSTEMS  DIVISION 
ATTN:  B  CUMMINGS 
ATTN:  C BRADFORD 

BDM  CORP 

ATTN:  J  LEECH 

BERKELEY  RSCH  ASSOCIATES,  INC 
ATTN:  S BRECHT 

BOEING  AEROSPACE  COMPANY 
ATTN:  NGERONTAKIS 

BOEING  TECHNICAL  &  MANAGEMENT  SVCS,  INC 
ATTN.  GHALL 

C.  L.  CONSULTING  SERVICES 
ATTN:  FFEER 

CALIF  RESEARCH  &  TECHNOLOGY,  INC 
ATTN:  M  ROSENBLATT 
ATTN:  RGAJ 
ATTN:  S KRUEGER 

CALSPAN  CORP 

ATTN:  R  MAMBRETTI 
ATTN:  R  MISSERT 

CARNEGIE  CORPORATION  OF  NEW  YORK 
ATTN:  D  ARSENIAN 

CARPENTER  RESEARCH  CORP 
ATTN:  HJ  CARPENTER 

CASSIDY  ANO  ASSOCIATES 
ATTN:  J JACOBSON 

CHARLES  SCAWTHORN 

ATTN.  C  SCAWTHORN 


CHARLES  STARK  DRAPER  LAB,  INC 
ATTN:  ATETEWSKI 


K 


COLORADO  STATE  UNIVERSITY 
ATTN:  D KRUEGER 
ATTN:  W  COTTON 

COMPUTER  SCIENCES  CORP 
ATTN:  G  CABLE 

DARTSIDE  CONSULTING 

ATTN:  A  FORESTER 

DEFENSE  ELECTRONICS  GROUP 
ATTN:  H  BURNSWORTH 
ATTN:  J  BELING 

DELTA  RESEARCH,  INC 

ATTN:  L  WEINER 
ATTN:  MRADKE 

DYNAMICS  TECHNOLOGY,  INC 
ATTN.  DHOVE 

ENW  INTERNATIONAL,  LTD 
ATTN:  J  CANE 

EOS  TECHNOLOGIES,  INC 
ATTN.  B GABBARD 
ATTN:  N JENSEN 
ATTN:  W  LELEVIER 

FACTORY  MUTUAL  RESEARCH  CORP 
ATTN:  M  A  DELICHATSIOS 
ATTN:  R  FRIEDMAN 

FEDERATION  OF  AMERICAN  SCIENTISTS 
ATTN:  J STONE 

GENERAL  ELECTRIC  CO 

ATTN:  RE  SCHMIDT 

GENERAL  RESEARCH  CORP 
ATTN:  B  BENNETT 
ATTN:  JBALTES 

HAROLD  ROSENBAUM  ASSOCIATES,  INC 
ATTN:  G  WEBER 

HORIZONS  TECHNOLOGY  INC 
ATTN:  A  EDWARDS 
ATTN:  J  AMBROSE 

HORIZONS  TECHNOLOGY,  INC 
ATTN.  RWLOWEN 
ATTN:  WTKREISS 

HUGHES  AIRCRAFT  CO 
ATTN:  EDIVITA 

IIT  RESEARCH  INSTITUTE 

ATTN:  H  NAPADENSKY 


Dist-4 


DNA-TR-87-176  (DL  CONTINUED) 


INFORMATION  SCIENCE,  INC 
ATTN:  W  DUDZIAK 

INSTITUTE  FOR  DEFENSE  ANALYSES 
ATTN:  C  CHANDLER 
ATTN:  E BAUER 
ATTN:  FALBINI 
ATTN:  L  SCHMIDT 

JOHNS  HOPKINS  UNIVERSITY 
ATTN:  M  LENEVSKY 
ATTN:  R  FRISTROM 
ATTN:  WBERL 

KAMAN  SCIENCES  CORP 
ATTN:  J  RUSH 
ATTN:  J SCRUGGS 

KAMAN  SCIENCES  CORP 
ATTN:  P  GRIFFIN 
ATTN:  P  TRACY 

KAMAN  SCIENCES  CORP 
ATTN:  E CONRAD 

KAMAN  SCIENCES  CORPORATION 
ATTN:  D ANDERSON 
ATTN.  DASIAC 

KAMAN  TEMPO 

ATTN:  B  GAMBILL 
ATTN:  D  FOXWELL 
ATTN:  DASIAC 
ATTN:  E  MARTIN 
ATTN:  R RUTHERFORD 
ATTN:  R  YOUNG 
ATTN:  SFIFER 
ATTN:  W  KNAPP 

LOCKHEED  MISSILES  &  SPACE  CO,  INC 
ATTN:  J  HENLEY 
ATTN:  J  PEREZ 

LOCKHEED  MISSILES  &  SPACE  CO.  INC 
ATTN:  P DOLAN 
ATTN:  W  MORAN 

M  I  T  LINCOLN  LAB 

ATTN:  S  WEINER 

MARTIN  MARIETTA  DENVER  AEROSPACE 
ATTN:  D  HAMPTON 

MAXIM  TECHNOLOGIES,  INC 
ATTN:  J MARSHALL 

MCDONNELL  DOUGLAS  CORP 
ATTN:  TCRANOR 
ATTN:  TTRANER 


MCDONNELL  DOUGLAS  CORP 
ATTN:  A  MONA 
ATTN:  F  SAGE 
ATTN:  G  BATUREVICH 
ATTN:  J GROSSMAN 
ATTN:  R  HALPRIN 
ATTN:  S  JAEGER 
ATTN:  WYUCKER 

MERIDIAN  CORP 

ATTN:  E  DANIELS 
ATTN:  F  BAITMAN 

MIDWEST  RESEARCH  INSTITUTE 
ATTN:  JS  KINSEY 

MISSION  RESEARCH  CORP 

ATTN:  R  ARMSTRONG 

MISSION  RESEARCH  CORP 
ATTN:  CLONGMIRE 
ATTN:  D ARCHER 
ATTN:  DKNEPP 
ATTN:  DSOWLE 
ATTN:  FFAJEN 
ATTN:  J  BALL 
ATTN:  KRCOSNER 
ATTN:  RBIGONI 
ATTN:  R  GOLDFLAM 
ATTN:  R  HENDRICK 
ATTN:  TOLD 
ATTN:  W  WHITE 

MITRE  CORPORATION 

ATTN:  J  SAWYER 

MODELING  SYSTEM,  INC 
ATTN:  G  BERLIN 

MRJ  INC 

ATTN:  D  FREIWALD 

NATIONAL  ADVISORY  COMMITTEE 
ATTN:  JALMAZAN 
ATTN:  J  BISHOP 

NATIONAL  INST.  FOR  PUBLIC  POLICY 
ATTN:  K  PAYNE 

NICHOLS  RESEARCH  CORP,  INC 
ATTN:  H  SMITH 
ATTN:  J  SMITH 
ATTN:  M  FRASER 
ATTN:  RBYRN 

NORTHROP  SERVICES  INC 
ATTN.  T  OVERTON 

NOTRE  DAME  DU  LAC,  UNIV  OF 
ATTN.  T  J  MASON 


>'f  k  ndi'T«i 


DNA-TR-87-176  (OL  CONTINUED) 

ORLANDO  TECHNOLOGY  INC 
ATTN:  R  SZCZEPANSKI 

PACIFIC-SIERRA  RESEARCH  CORP 
ATTN.  GANNO 

ATTN:  H  BRODE,  CHAIRMAN  SAGE 
2  CYS  ATTN:  K  HEIKES 
2  CYS  ATTN  L  RANSOHOFF 
ATTN:  MDORE 
2  CYS  ATTN:  R  SMALL 

PHOTOMETRICS,  INC 

ATTN:  I  L  KOFSKY 

PHOTON  RESEARCH  ASSOCIATES 
ATTN:  JMYER 

PHYSICAL  RESEARCH  CORP 
ATTN:  ACECERE 

PHYSICAL  RESEARCH  INC 
ATTN:  HFITZ 

PHYSICAL  RESEARCH  INC 
ATTN:  D  MATUSKA 

PHYSICAL  RESEARCH  INC 
ATTN:  J  WANG 
ATTN:  WSHIH 

PHYSICAL  RESEARCH,  INC 

ATTN.  D  WESTPHAL 
ATTN:  DWHITENER 
ATTN:  H  WHEELER 
ATTN:  R  BUFF 
ATTN:  R  DELIBERIS 
ATTN:  T  STEPHENS 
ATTN:  WC  BLACKWELL 

PHYSICAL  RESEARCH,  INC 
ATTN:  G HARNEY 
ATTN:  J  DEVORE 
ATTN:  J  THOMPSON 
ATTN:  R  STOECKLY 
ATTN:  WSCHLUETER 

PHYSICAL  RESEARCH,  INC 
ATTN:  H  SUGIUCHI 

POLYTECHNIC  OF  NEW  YORK 
ATTN:  BJ  BULKIN 
ATTN:  GTESORO 

PRINCETON  UNIVERSITY 

ATTN.  J  MAHLMAN 

QUADRI  CORP 

ATTN:  HBURNSWORTH 


R  &  D  ASSOCIATES 

ATTN:  AKUHL 
ATTN:  D  HOLLIDAY 
ATTN:  F  GILMORE 
ATTN:  G  JONES 
ATTN:  J SANBORN 
ATTN:  RTURCO 

R  &  D  ASSOCIATES 

ATTN:  B  YOON 

R  J  EDWARDS  INC 

ATTN:  R  SEITZ 

RADIATION  RESEARCH  ASSOCIATES,  INC 
ATTN:  B  CAMPBELL 
ATTN:  M  WELLS 

RAND  CORP 

ATTN:  GL  DONOHUE 
ATTN.  P  DAVIS 
ATTN:  P  ROMERO 

RAND  CORP 

ATTN:  B  BENNETT 
ATTN:  J  GERTLER 

ROCKWELL  INTERNATIONAL  CORP 
ATTN:  J  KELLEY 

S-CUBED 

ATTN.  B  FREEMAN 
ATTN:  KDPYATT,  JR 
ATTN:  R  LAFRENZ 

S-CUBED 

ATTN:  C NEEDHAM 
ATTN:  SHIKIDA 
ATTN:  T CARNEY 

SCIENCE  APPLICATIONS  INC 
ATTN:  REDELMAN 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  CHILL 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN.  B  MORTON 
ATTN:  B  SCOTT 
ATTN:  D SACHS 
ATTN:  GT  PHILLIPS 
ATTN:  JBENGSTOM 
ATTN:  M  DRAKE 
ATTN:  M  MCKAY 
ATTN:  D  HAMLIN 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  D BACON 
ATTN:  DRLGOURE 
ATTN:  F  GIESSLER 
ATTN:  J COCKAYNE 


Dist-6 


m 


IS* 


SCIENCE  APPLICATIONS  INTL  CORP 
ATTN.  J  SONTOWSKI 


SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  T  HARRIS 

SCIENTIFIC  RESEARCH  ASSOC,  INC 
ATTN:  B  WEINBERG 

SCIENTIFIC  SERVICES,  INC 
ATTN:  C  WILTON 

SRI  INTERNATIONAL 

ATTN:  CWITHAM 
ATTN:  D  GOLDEN 
ATTN:  D  MACDONALD 
ATTN:  D ROBERTS 
ATTN:  EUTHE 
ATTN:  GABRAHAMSON 
ATTN.  J  BACKOVSKY 
ATTN:  W  CHESNUT 

SRI  INTERNATIONAL 

ATTN:  RBRAMHALL 
ATTN:  R  WOOLFOLK 
ATTN:  W  VAIL 

STAN  MARTIN  ASSOCIATES 
ATTN:  SB  MARTIN 
ATTN:  S  MARTIN 

STANTON  CONSULTING 

ATTN:  M STANTON 

SWETL.  INC 

ATTN:  T  PALMER 
ATTN:  TY  PALMER 

SYSTEM  PLANNING  CORP 
ATTN:  B  GARRETT 
ATTN:  C  FELDBAUM 
ATTN:  J  SCOURAS 
ATTN:  MBIENVENU 
ATTN:  RSCHEERBAUM 

SYSTEMS  AND  APPLIED  SCIENCES  CORP 
ATTN:  M  KAPLAN 

TECHNOLOGY  INTERNATIONAL  CORP 
ATTN:  W  BOQUIST 

TELEDYNE  BROWN  ENGINEERING 
ATTN:  AORTELL 
ATTN:  F  LEOPARD 
ATTN:  J  FORD 


"  . 
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TELEDYNE  BROWN  ENGINEERING 
ATTN:  DGUICE 


TEXAS  ENGR  EXPERIMENT  STATION 
ATTN:  W  H  MARLOW 


TOYON  RESEARCH  CORP 


ATTN 

ATTN 

ATTN 


C  TRUAX 
J  GARBARINO 
J  ISE 


TRW  ELECTRONICS  &  DEFENSE 

ATTN: 

M  HAAS 

TRW  INC 

ATTN: 

F  FENDELL 

ATTN: 

G  KIRCHNER 

ATTN: 

H  CROWDER 

ATTN: 

J  FEDELE 

ATTN: 

M  BRONSTEIN 

ATTN: 

R  BACHARACH 

ATTN: 

S  FINK 

ATTN. 

T  NGUYEN 

VISIDYNE,  INC 

ATTN: 

H  SMITH 

ATTN: 

J  CARPENTER 

WASHINGTON, 

UNIVERSITY  OF 

ATTN:  J  I  KATZ 


FOREIGN 


AERE  ENVIRONMENTAL  AND  MEDICAL  SC 
ATTN:  SPENKETT 


ATOMIC  WEAPONS  RESEARCH  ESTABLISHMENT 
ATTN:  PF  A  RICHARDS 


ATOMIC  WEAPONS  RESEARCH  ESTABLISHMENT 
ATTN:  DL  JONES 
ATTN:  D  M  MOODY 


AUSTRALIA  EMBASSY 

ATTN:  DR  LOUGH 
ATTN:  MAJ  GEN  H  J  COATES 
ATTN:  P  PROSSER 


BRITISH  DEFENCE  STAFF 
ATTN:  C  FENWICK 
ATTN: 

ATTN: 

ATTN: 

ATTN: 


J  CRANIDGE 
J  EDMONDS 
M  NORTON 
P  WEST 


CANADIAN  FORESTRY  SERVICE 
ATTN:  B  STOCKS 
ATTN:  TLYNHAM 


CSIRO 


ATTN:  IGALBALLY 


$ 


a 


$ 


'4 


v 


53 


Dist-7 


DNA-TR87-176  (DL  CONTINUED) 


CSIRO:  ATMOSPHERIC  RESEARCH 
ATTN:  APITTOCK 


HARVARD  UNIVERSITY 

ATTN:  G  CARRIER 


EMBASSY  OF  BELGIUM 

ATTN.  LARNOULD 

ISRAEL  EMBASSY 

ATTN.  N  BELKIND 

MAX-PLANCK  INSTITUTE  FOR  CHEMISTRY 
ATTN:  P  J  CRUTZEN 


MINISTRY  OF  DEFENCE 
ATTN:  R  RIDLEY 


HARVARD  UNIVERSITY 

ATTN:  DEARDLEY 


i 

IOWA,  UNIVERSITY  OF 

ATTN.  HISTORY  DEPT/S  PYNE 


MARYLAND  UNIVERSITY  OF 

ATTN:  A  ROBOCK  DEPT  METEOROLOGY 
ATTN:  A  VOGELMANN  DEPT  METEOROLOGY 
ATTN:  R  ELLINGSON  DEPT  METEOROLOGY 


NATIONAL  DEFENCE  HEADQUARTERS 
ATTN:  HAROBITALLE 

TRINITY  COLLEGE 

ATTN:  F  HARE 

DIRECTORY  OF  OTHER 


MIAMI  LIBRARY  UNIVERSITY  OF 
ATTN:  C CONVEY 

MIAMI  UNIV  LIBRARY 

ATTN.  J  PROSPERO  ATMOS  SC 

NEW  YORK  STATE  UNIVERSITY  OF 
ATTN:  R  CESS 


ATMOS.  SCIENCES 

ATTN:  GSISCOE 


OAK  RIDGE  ASSOCIATED  UNIVERSITIES 
ATTN:  C  WHITTLE 


BROWN  UNIVERSITY 

ATTN:  RK  MATTHEWS 

BUCKNELL  UNIVERSITY 

ATTN:  0  ANDERSON 


PENNSYLVANIA  STATE  UNIVERSITY 
ATTN:  DWESTPHAL 

SOUTH  DAKOTA  SCH  OF  MINES  &  TECH  LIB 
ATTN:  H  ORVILLE 


CALIFORNIA,  UNIVERSITY 

ATTN:  R  WILLIAMSON 


TENNESSEE,  UNIVERSITY  OF 
ATTN:  K  FOX 


CALIFORNIA,  UNIVERSITY  OF 

ATTN.  L  BADASH/DEPT  OF  HIST 


UNIVERSITY  OF  SOUTH  FLORIDA 
ATTN:  SYING 


COLORADO,  UNIVERSITY  LIBRARIES 
ATTN:  JBIRKS 
ATTN:  R  SCHNELL 

DREXEL  UNUVERSITY 

ATTN:  J  FRIEND 

DUKE  UNIVERSITY 

ATTN:  F  DELUCIA 

GEORGE  MASON  UNIVERSITY 
ATTN  PROFS  SINGER 
ATTN:  R  EHRLICH 

GEORGE  WASHINGTON  UNIVERSITY 
ATTN:  R  GOULARD 

GEORGIA  INST  OF  TECH 

ATTN:  E  PATTERSON 

HARVARD  COLLEGE  LIBRARY 
ATTN:  W  PRESS 


UNIVERSITY  OF  WASHINGTON 
ATTN:  CLEOVY 
ATTN:  LRAOKE 
ATTN:  P  HOBBS 

VIRGINIA  POLYTECHNIC  INST  LIB 
ATTN:  M  NADLER 

WASHINGTON  STATE  UNIVERSITY 
ATTN:  DR  A  CLARK 

WISCONSIN  UNIVERSITY  OF 
ATTN:  P  WANG 


I 

I 
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